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A similarity hypothesis for the 
two-point correlation tensor in a 
temporally evolving plane wake 


By D. W. Ewing 1 , W. K. George 2 , R. D. Moser 3 , AND M. M. Rogers 4 


1. Motivation and objectives 

It has long been known that the equations that govern the evolution of the single- 
point moments, such as the mean velocity or the turbulent Reynolds stresses, admit 
similarity solutions for many of the basic shear flows ( e.g . George 1989 or Tennekes 
and Lumley 1972). In this approach, it is argued that the flow evolves such that 
all of the terms in the governing equations make the the same relative contribution 
so the flow reaches an equilibrium or similarity state. In many cases the initial 
conditions of the flow are inconsistent with the hypothesized similarity solutions so 
these solutions are, at most, an approximation of the flows asymptotic state. How- 
ever, the agreement between the predictions of the theory and experimental data 
(e.g., Wygnanski et a l. 1989) suggests that flows do approach such an asymptotic 
state. 

Traditionally (e.g., Tennekes and Lumley 1972), it was argued that this asymp- 
totic state is universal for all flows of a particular type (e.g., all plane wakes). 
However, this argument is not consistent with the measurements in the far field of 
plane wakes reported by Wygnanski et al. (1989), who found that the similarity 
profiles of the normal Reynolds stresses (particularly the streamwise component) 
differed for different types of wake generators. George (1989) argued that this oc 
curred because the governing equation for the similarity profile of the turbulent 
kinetic energy contained a constant that depended on the growth rate of the flow. 
Consequently, George (1989) concluded that the similarity profiles for the normal 
stresses would not be universal unless the growth rates of all wakes were the same, 
which was not the case for the wakes studied by Wygnanski et al. (1986). (A simi- 
lar conclusion was reached from a later analysis of the governing equations for the 
individual normal stresses; e.g., George 1994 or Ewing and George 1994.) George 
(1989) attributed the differences in the growth rates of the wakes (and hence the 
asymptotic states) to differences in the coherent structures produced by the bodies 
generating the wake. 
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It is now widely recognized that coherent structures play an integral role in many 
of the processes that determine the growth rate of a turbulent flow, including the 
entrainment of irrotational flow and the gross mixing of the fluid across the layer. 
A question that has not yet been resolved is whether the structures in the flow 
themselves (or more properly the coherent structures and the probability density 
functions describing their occurrences) approach an ‘equilibrium’ or similarity state 
when the single-point moments measured in the flow agree with the predictions 
of the similarity hypothesis. That is, can all of the statistical measures of the 
structures (the single-point moments being the simplest) at different points in a 
flow’s evolution be related by a similarity transformation. 

The primary objective of this research is to examine whether the governing equa- 
tions for more complex statistical measures of the structures in the temporally 
evolving plane wake admit similarity solutions. (In the initial stage of this research 
it was established that the governing equations for the mean momentum and the 
Reynolds stresses admit similarity solutions, v. Moser et al. 1995). The two-point 
velocity correlation tensor was chosen because it contains more information about 
the turbulent structure than the single-point moments and it is often used in at- 
tempts to educe coherent structures from the flow (e.y., Grant 1958 or Payne and 
Lumley 1967). There have been few previous attempts to demonstrate that the 
governing equations for these two-point correlations admit similarity solutions in 
non-homogeneous flows. Ewing and George (1994) and Ewing (1995) demonstrated 
that the governing equations for the two-point velocity correlation tensor in the far 
fields of the spatially evolving axisymmetric and plane jets admit similarity solu- 
tions. However, in both of these cases the predictions of the similarity hypotheses 
were not tested using experimental data, so it was not determined if these similarity 
solutions are accurate descriptions of real shear flows. This question is addressed 
in this research using data from Direct Numerical Simulations of the temporally 
evolving wake computed by Moser and Rogers (1994). 


2. Accomplishments 

2.1 Theoretical analysis 

In the temporally evolving wake, shown in Fig. 1, a momentum deficit spreads 
in one non-homogeneous direction, X 2 , as the flow evolves in time. The other 
two directions, including the direction of the mean flow, x u are homogeneous. 
Consequently, the governing equations for the two-point velocity correlations in 
this flow can be written as 
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Figure 1. Geometry of the temporally evolving wake 
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where C = ii — x\, 7 = x 3 - x ' 3 , U\ is the mean velocity in the 27 -direction, and u, 
are the fluctuating velocity components in the 27 directions. In these equations, the 
primed variables are evaluated at one arbitrary point in space while the unprimed 
variables are evaluated at a second arbitrary point. The two-point velocity corre- 
lations must also satisfy initial conditions and a number of boundary conditions; 


r 2 ,* 2 ,* , 7— ' >±00 


Uiu'j = 0 , 


assuming the turbulence in the free stream is negligible. 

Following the approach outlined by George (1989), it is hypothesized that the 

governing equations for the two-point correlation tensor admit similarity solutions 
given by 

«.-(*!, *2, *3,0 = Q iJ (t)qiAt, v, v', C> *), (3a) 

* x 1^2,X3,0^'„x',l',t) = PfWp'jtM 1 , C,*), (3b) 

P'(*i,*2,*3,<W*i,*2, *3,0 = V, C, *), (3c) 

UkUiu) = T^' 3 {t)U\ i} (^ t), V, c, *), (3d) 

u ' u k u 'j = T2 k *(t)ttl,kj(Zi r h T )\C,*)’> (3e) 
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and the * in Eqs. 3a -3e indicate that the solutions may depend on the source 
conditions of the flow. At this point the length scales 6 X and S 3 are arbitrary. 
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The allowable choices for these scales are determined by examining the equations 
of motion. The length scale used to scale the x 2 is equal to the scale used in the 
single-point analysis, 6 (t) oc (t — t 0 ) 1 ' 2 (Moser et al. 1995), since the similarity 
solutions for the single- and two-point correlations must be consistent in the limit 
when the separation distance between the two points is zero. Here, t a is the location 
of the virtual origin of the wake. 

Substituting these hypothesized solutions for the two point correlations and the 
similarity solution for the mean velocity given by (Moser et al. 1995) 

Uloo - Ui(x 2 ,t) = Us{t)f{r]) (4) 


into Eq. 1 yields 



where the time-dependent portion of each term is contained in square brackets. 
It is possible to remove the time dependence from these equations if the time- 
dependent portion of the all terms in each equation are proportional, leaving a 
set of equations for the similarity solutions in terms of the similarity variables and 
the constants whose values depend on the ratios of the time-dependent terms. In 
this case, the similarity solutions again represent an ‘equilibrium’ solution for the 
governing equations since all of the terms make the same relative contribution as 
the flow evolves. 

It is straightforward to demonstrate that the time dependent portions of the four 
viscous terms in Eq. 5 are only proportional if 

oc 6 oc 63 oc (t — to) 1 / 2 . 

This implies that all three length scales must have the same virtual origin. 


( 6 ) 



Similarity solutions for multi-point correlation 


167 


Further, the time-dependent portions of the convective and the unsteady terms 
in Eq. 5 are only proportional to the time-dependent portions of these viscous terms 
if 


'Q^dS' 
6 dt 


oc 




(7a) 


and 



(7b) 


These constraints are satisfied if the three time scales 6/U a (a characteristic time 
scale of the mean strain rate), S 2 /u (a viscous time scale), and (1/6 d6/dt ) _1 (a 
charteristic time scale of the spreading of the flow) are proportional. Thus, the ratio 
of these time scales are constant; i.e., 
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and 
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which are the constraints derived in the analysis of the single-point equations. 

The constraint that the Reynolds number is a constant is a gratifying result since 
a single length scale was used to define each of the similarity coordinates in Eq. 3f. 
Thus, it was implicitly assumed that all of the physically relevant length scales in the 
flow grow in proportion as the flow evolves; an assumption that is generally thought 
to be valid in a constant-Reynolds-number flow (Batchelor 1953). The functional 
form of the solutions will, of course, depend on the ratios of the physically relevant 
length scales (or the Reynolds number). 

The constant /3 that appears in the equations is analogous to the constant that 
George (1989) found in the governing equation for the turbulent kinetic energy in 
the spatially evolving wake. It is very likely that the value of this ratio depends 
on the coherent structures in the flow, and it is through this constant that their 
influence is incorporated into the equations for the similarity solutions. 

The time-dependent portions of the rest of the terms in Eq. 5 are proportional if 


Pi <xU.Q i k ,Pi j oc U,Q k >, 


and 

T k, ' j cc T 2 ' k] oc U s Q' j . (9b) 

These expressions must be satisfied for any choice of k. The choice for Q'' 1 is not 
uniquely determined from Eq. 5; however, Q' 1 is not arbitrary since the similarity 
solution for the two-point velocity correlation tensor must be consistent with the 
similarity solution derived for the single-point moments (Moser et al. 1995). Thus, 
it follows that 


Q i J oc ui 


( 10 ) 
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Of course, these constants of proportionality may be functions of 0 or Res- 

Thus, the governing equations for the two-point velocity correlation tensor in the 
temporally evolving wake admit similarity solutions, of which the similarity solution 
for the single-point Reynolds stresses are a special case. It is straightforward to 
demonstrate that the governing equations for the similarity solutions include the 
constants Re s and 0. Thus, in general, the similarity solutions are functions of these 
two ratios whose values depend on the source conditions of the flow. In many flows, 
the initial conditions are not consistent with the hypothesized similarity solution so 
these solutions axe, at most, an approximation of the asymptotic state of the flow 
(as with the single-point similarity solutions). 

2.2 Implications of the similarity hypothesis 

When a similarity solution exists for the two-point velocity tensor, other statis- 
tical measures that can be determined directly from this two-point correlation also 
have similarity solutions. In many cases, these results provide useful predictions to 
compare with data in order to test the similarity hypothesis. 

For example, when a similarity solution exists for the two-point velocity correla- 
tion tensor, the one-dimensional spectra in the x \ -direction given by 

1 f°° 

Flj(ki,x 2 ,x' 2 ,7) = — j Rij(<;,x 2,x 2 ,7) e 1 


(where R t J = Uiu'j) can be written as 

Fij(kuX2,X2,7) = F^(fc (11a) 

where the similarity solution for the one- dimensional spectra is given by 

= 2^ y_ oo 9.,iU,r / ,r / ',C)e- , * 1 «de (He) 


and k\ is the given by 

fcj = kiS. (Hb) 

A similar relationship can be derived for the one- dimensional spectra in the re- 
direction. Thus, when a similarity solution exists for the two-point velocity correla- 
tion tensor, the one-dimensional spectra of the field maintain the same shape while 
continuously shifting downward and to lower wavenumbers as the flow spreads. 

It is also straightforward to demonstrate that many of the classically defined 
turbulent length scales are proportional to the similarity length scale. For example, 
the Taylor microscales given by 
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can be written as (Ewing 1995) 


X i,f) = (13) 

It is also straightforward to show that the integral length scales and the Kolmogorov 
length scale are proportion to the similarity length scale (Ewing 1995). 

A number of other interesting implications of the two-point similarity hypothe- 
sis, such as similarity solutions for the moments related to the two-point velocity- 
gradient correlation tensor (e.g., the two-point vorticity correlation or the dissipa- 
tion of the turbulent kinetic energy in the flow), are discussed in Ewing (1995) and 
Ewing et ai (1995). 


2.3 Comparison with simulation data 


Although, the previous analysis demonstrated, that the hypothesized similarity 
solutions are consistent with the equations of motion, this does not ensure that these 
solutions are an accurate description of real flows. In this case the predictions of 
the similarity hypothesis for the two-point velocity correlation tensor are compared 
to the data from two Direct Numerical Simulations of the temporally evolving wake 
computed by Moser and Rogers (1994). These simulations were carried out in finite 
boxes with periodic boundary conditions in the homogeneous directions in contrast 
to the theory, which was developed for a wake in an infinite domain. Consequently, 
the theory cannot exactly approximate the flow in these simulations. However, it is 
generally argued that the evolution of the scales of motion that are ‘much 1 smaller 
than the dimension of the boxes in the periodic directions should be similar to the 
evolution of those in an infinite wake. Thus, a comparison between the predictions 
of the similarity hypothesis and the data from the simulations is essentially a test 
of both the similarity hypothesis and this idea. Agreement between the predictions 
of the theory and the data lends support to both. (Experimentalists experience the 
very same problem when data from finite experimental rigs, such as wind tunnels, 
are used to test hypotheses developed for infinite flow; e.o., George and Gibson 
1992). 


The initial conditions for the wake simulations were generated using two realiza- 
tions from a turbulent boundary layer simulation yielding a wake with a Reynolds 
number, given by 




S?JVi - U loo )dx 2 


(14) 


of 2000. In the first (unforced) wake simulation the base initial conditions were 
used to initiate the flow while in the second (forced) wake simulation the uj and 
u 2 velocity components of the spanwise two-dimensional modes were amplified by 
a factor of 5. These simulations are computed in periodic boxes of length 50^ in 
the xi -direction and 12.5 6 d in the x 3 -direction, where 6 d is the initial displacement 
thickness of the wake. Moser and Rogers (1994) found that the data from both 
simulations were approximately in agreement with the predictions of the similarity 
hypothesis for a period of the flow s evolution. For example, it is evident from 
Figs. 2a and 2b that the Reynolds number of both flows is approximately constant 
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FIGURE 2. Evolution of the Reynolds number and similarity length scale: : 

unforced wake simulation; : forced wake simulation; o : location of points 

used to examine the two-point similarity hypothesis. 


and the square of the similarity length scales in the flows grow approximately lin 
early, as predicted in the similarity analysis, for a significant time period. The 
two-point similarity hypothesis is tested using data from four points in this range, 
indicated by the diamonds in Figs. 2a and 2b. 

The unsealed one- dimensional spectra F/j raid at the centerline of the two 
wakes (spatially averaged in the X 3 -direction) are illustrated in Figs. 3a and 3b. 
The spectra exhibit peaks in the low-wavenumber region that are inconsistent with 
the similarity hypothesis since they occur at the same wavenumber in physical vari- 
ables. However, it is not anticipated that the similarity hypothesis should collapse 
the spectra in this region because it is likely that these motions are affected by 
the periodic boundary conditions or, conversely, the coarse discretization of wave 
space at these scales. In contrast, the spectra in the high- wavenumber region shift 
downward and to the left as the flow evolves, in agreement with the predictions of 
the similarity hypothesis. The amplitude of the spectra in this high- wavenumber 
region differ by a factor of approximately 3 — 4 so they should provide a good test 
of the similarity hypothesis. 

The scaled one- dimensional spectra F ^ at the centerline and the half-deficit 
point, 77 = 0.5, in the unforced wake are illustrated in Fig. 4, while the scaled 
one-dimensional spectra Fk at the same points in the forced wake simulation are 
illustrated in Fig. 5. In both of these figures the data from the half-deficit point 
in the wake are shifted up by a order of magnitude. Overall there is excellent col- 
lapse of the data for the region where fci > ISAfci, indicating that the statistical 
measures of all but the largest motions evolve as predicted by the similarity hypoth- 
esis. There is some discrepancy between the predictions and the data at the largest 
wavenumbers; however, these variations occur because the effective resolution of the 
simulations varies as the flow evolves (Ewing 1995). The one- dimensional spectra 
of the correlations in the ^-direction (Ewing 1995) exhibit better collapse than the 
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FIGURE 3. Unsealed one-dimensional spectra in (a) the unforced wake, (b) the 

forced wake: : F/j, : F\ 2 . In (a) 6/6 d is equal to: 3.494 o , 4.008 □ , 

4.446 o , and 4.884 A . In (b) S/6 d is equal to: 5.306 o , 6.037 □ , 6.636 o , and 7.235 



FIGURE 4. Scaled one-dimensional spectra, F / 3 in the unforced wake at rj = 

0.0 and r/ = 0.5. 6/6 d is equal to; 3.494 o , 4.008 □ , 4.446 o , and 4.884 a . 

spectra in the x \ -direction (partially because these spectra are spatially averaged 
in the x\ -direction, which is 4 times longer than the :r 3 -direction used to average 
F 1 ). 

The large-scale motions make a contribution to the two-point velocity correlation 
tensor for all separation distances. Consequently, the structure functions in the 
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FIGURE 5. Scaled one-dimensional spectra, F 22 in the forced wake at 77 = 0.0 

and 77 = 0.5. 6/6d is equal to; 5.306 o , 6.037 a , 6.636 o , and 7.235 ^ . 

r 2 -direction, given by 


(u a - u' Q ) 2 — u 2 a + u' 2 - 2R a , Q (x2,X2) (15) 

are used to compare the prediction of the similarity hypothesis and the data, since 
they are more measures of the motions whose sizes are on the order of the separation 
distance between the points. The structure functions for a = 1 centered around the 
centerline and the half-deficit point in the unforced wake, seeded with appropriate 
similarity variables, axe shown in Figs. 6a and 6b. In both cases the profiles from 
the four different times collapse for small and intermediate separation distances, 
suggesting again that the statistical measures of all but the largest scales of motion 
in the flow are evolving as predicted by the similarity hypothesis. The structure 
functions from the forced wake simulation also collapse for small and intermediate 
separation distances when they are scaled using similarity variables (Ewing 1995). 

2-4 Scalar fields 

The similarity analysis was also extended to the governing equations for a passive 
scalar field. It was demonstrated that these equations admit similarity solutions for 
both a scalar field with a mean deficit in the wake region and a scalar field where 
the mean value in the two free streams differ. The predictions of the similarity 
hypothesis were compared to data from the simulations, which were computed with 
this second type of scalar field. However, the variations of the two-point scalar 
correlation in the simulations were much smaller than the variations of the two- 
point velocity correlation tensor (as predicted by the similarity hypothesis). As 
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FIGURE 6. Scaled structure functions about (a) Tf = 0 (b) rj = 0.5. 6/6d is equal 
to: 3.494 o , 4.008 □ , 4.446 o , and 4.884 ^ . 

a result, the data from the period where the simulations were approximately self- 
similar did not vary sufficiently to rigorously test the predictions of the similarity 
hypothesis for the scalar field. Further details of this comparison can be found in 
Ewing (1995). 

3. Conclusions 

The analysis demonstrated that the governing equations for the two-point veloc- 
ity correlation tensor in the temporally evolving wake admit similarity solutions, 
which include the similarity solutions for the single- point moment as a special case. 
The resulting equations for the similarity solutions include two constants, f3 and 
Res , that are ratios of three characteristic time scales of processes in the flow: a 
viscous time scale, a time scale characteristic of the spread rate of the flow, and a 
characteristic time scale of the mean strain rate. The values of these ratios depend 
on the initial conditions of the flow and axe most likely measures of the coher- 
ent structures in the initial conditions. The occurrences of these constants in the 
governing equations for the similarity solutions indicates that these solutions, in 
general, will only be the same for two flows if these two constants are equal (and 
hence the coherent structures in the flows axe related). 

The compaxisons between the predictions of the similarity hypothesis and the data 
presented here and elsewhere (Ewing 1995) indicate that the similarity solutions for 
the two-point correlation tensors provide a good approximation of the measures of 
those motions that are not significantly affected by the boundaxy conditions caused 
by the finite extent of real flows. Thus, the two-point similarity hypothesis provides 
a useful tool for both numerical and physical experimentalist that can be used to 
examine how the finite extent of real flows affect the evolution of the different scales 
of motion in the flow. The similarity analysis of the governing equations for the 
multi-point correlations can be extended to a wide range of spatially and temporally 
evolving flows using the methodology outlined by Ewing (1995), so this technique 
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can be used to examine the effect of finite boundaries on the evolution of a number 
of different flows. 


REFERENCES 

Ewing, D. & George, W. K. 1994 Applications of a similarity hypothesis to the 
proper orthogonal decomposition for spatially evolving flows. In Proc . of the 
Int . Symp. on Turb. Heat and Mass Trans ., Lisbon, Portugal. 

Ewing, D. 1995 On multi-point similarity solutions in turbulent free-shear flows. 
Ph.D. Dissertation, State University of New York at Buffalo. 

Ewing, D., George, W. K., Rogers, M. M.& Moser, R. D. A similarity 
hypothesis for the two-point correlations in a temporally evolving, plane wake 

J. Fluid Meek, in preparation. 

GEORGE, W. K. 1989 The self-similarity of turbulent flows and its relation to 
initial conditions and coherent structures. In Advances in Turbulence, edit. W. 

K. George and R. E. Arndt, Hemisphere Publishing. 

GEORGE, W. K. & Gibson, M. M. 1992 The self-preservation of homogeneous 
turbulence. Expt . in Fluids. 13, 229. 

GEORGE, W. K. 1994 Some new idea for similarity of turbulent shear flows. In 
Proc. of the Int. Symp. on Turb. Heat and Mass Trans., Lisbon, Portugal. 

GRANT, H. L. 1958 The large eddies of turbulent motion. J. Fluid Mech . 3, 149. 

MOSER, R. D. & Rogers M. M. 1994 Direct simulation of a self-similar plane 
wake. NASA Tech. Memo, TM 108815. 

Moser, R. D., Rogers, M. M., & Ewing, D. 1995 J. Fluid Mech. in preparation. 

PAYNE, F. R. &; LUMLEY, J. L. 1967 Large eddy structure of the turbulent wake 
behind a circular cylinder. Phys. Fluids Supp . 10, 194. 

Tennekes, H. & Lumley, J. L. 1972 A First Course in Turbulence. MIT Press. 

Wygnanski, I., Champagne, F., & Marasli, B. 1986 On the large-scale struc- 
tures in two-dimensional, small- deficit, turbulent wakes. J. Fluid Mech. 168, 
31. 



